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,S^ • I review the array of event generators which have been written 
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I review the array of event generators which have been written to provide simula- 
tions of high-energy e+e" reactions. 



1 Introduction 

In any simulation of an experimental analysis at an e~^e~ collider, we must 
begin with a sample of physics events to be analyzed. To produce these, we 
need an event generator. This program encodes our knowledge of Standard 
Model background processes and our expectations for signal processes. In 
this article, I will review the current variety of event generators available for 
simulations studies at future linear colliders. 

There are a number of goals that an event generator might be expected 
to fulfill. It should realistically represent possible signal processes and Stan- 
dard Model backgrounds. It should take care of the superposition of QCD and 
fragmentation effects onto electroweak cross sections. It should give high ac- 
curacy, for precision studies. And, it should have the flexibility to include new 
reactions of arbitrary and exotic form. These goals generally conflict with one 
another, or else are achieved only at the expense of a high level of complexity. 
So we needed different tools optimized for these various tasks. 

In this report, I review the event generators which address these issues 
that were presented at Sitges. The major simulation programs described in 
this article are listed in Table H. This table includes, for each program, a Web 
address where download information and documentation can be found. 

2 Workhorses 



Among event generators, first place_must be given to the general purpose 
programs PYTHIA El and HERWIGH-H Both programs were originally written 
to test ideas about QCD jet phenomena and hadronization. But both have 
now evolved into general-purpose codes incorporating all of the basic Standard 
Model processes in e~^e~ annihilation and a variety of nonstandard reactions. 
The most important aspect of PYTHIA and HERWIG is that they fully 
simulate QCD final state state effects. Given a system of two or more par- 
tons with large invariant mass, these programs generate a QCD parton shower 
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Table 1: Event generators for e+e linear collider physics 

PYTHIA www . thep . lu . se/~torb j orn/Pythia . html 

HERWIG hepwww . rl . ac . uk/theory/seymour/herwig/ 

CIRCE heplix. ikp.physik.tu-darmstadt .de/nlc/beam.html 

PHYSSIM www-jlc .kek. jp/subg/of f 1/physsim/ 

PANDORA www. slac . Stanford. edu/~mpeskin/LC/pajidora.htinl 

ISAJET quark.phy.bnl.gov/~paige/ 

SUSYGEN lyoinf o . in2p3. f r/susygen/susygen3 .html 

EXCALIBUR home . cern . ch/charlton/excalibur/excalibur . html 

KORALW hpjmiady.if j .edu.pl/programs/node9.html 

WPHACT www.to.infn. it/'ballestr 

KK home . cern. ch/~jadach/ 

GRACE www-sc.kek.jp/minami/ 

COMPHEP theory .npi .msu. su/~comphep/ 



and then simulate the hadronization of the final array of partons. The shower 
algorithm is not exact at higher orders in a^ but does generate an approxi- 
mately correct set of hard jets radiated from the original parton system. The 
hadronization step is carried out by different algorithms in the the two pro- 
grams, but, in both cases, the description of hadronization has been tuned 
to fit the e+e~ annihilation data. These features imply that QCD final-state 
interactions have been included in a way that extrapolates correctly to high 
energy. 

PYTHIA allows any parton-level generator to be included as a hard sub- 
process. The generator must specific the color routing in the final state and 
the order in which parton showers are to be generated. One caution with this 
approach is that order a^ corrections can raise or lower the overall normaliza- 
tion of the cross section. This effect cannot be included in the hadronization 
but must be accounted for externally. 

PYTHIA can run with a given initialization at a variety of e+e~ center 
of mass energies. This allows the program to be linked to ^ generator of 
initial-state e~ and e"*" energy distributions, such as CIRCE Q or PYBMSp 
to simulation the effect of beamstrahlung. Initial state polarization is not 
included in the current version. Final state spin correlations are included for 
some but not all processes; notably, spin correlations are included for the very 
important process e+e^ -^ W^W^ . 

PYTHIA and HERWIG both include generators for the process 77 — > 
hadrons, including hard, soft, and 'resolved' components. The third of these 



refers to processes in which partons in the photon undergo a hard scatter- 
ing. The relative magnitudes of these three components are not understood 
from theory, but it is important to understand this problem to compute the 
high-rate 'minijei' background in which a 77 collision produces a low-mass 
hadronic systemaQ New data on high energy 77 processes from LEP 2 and 
on 7p processes from HERA — which contain much of the same physics-should 
allow systematic tuning of these generators. 

3 Polarization 

To introduce the next sections of this report, I must digress on the subject 
of polarization. Polarization has a central role in the LC physics. On one 
hand, because the LC will operate in the energy region well above the Z" 
where it becomes manifest that left- and right-handed have completely different 
quantum numbers, all standard and non-standard cross sections will depend 
strongly on polarization. On the other hand, since it is difficult to measure 
polarization effects in the hadronic environment, polarization provides many 
new observables that cannot be studied at the LHC. Thus it is important 
that both initial- and final-state polarization be included properly in LC event 
generators. 

How can polarization be included in physics simulations? The traditional 
approach is to generalize cross section formulae to include polarization asym- 
metries. However, this rapidly becomes cumbersome. Generators that take 
polarization seriously typically work at the amplitude level. Even if one does 
not include polarization, it is useful to work with amplitudes in any complex 
Feynman diagram computation, since if there are N terms in the expression 
for the amplitude, there are iV^ in the expression for the cross section. Any 
method that makes use of amplitudes to compute the cross section can be 
structured so that the polarization-dependence is available for free. 

There are two common approaches for including polarization in the cross- 
section formulae used in event generators. The first approach is the helicity- 
amplitude paradigm. In this approach, one computes amplitudes for transitions 
between states of definite hclicityQ These amplitudes are then linked together 
to provide the complete amplitude for a process that turns the initial e+e^ 
state into the final decay products. Finally, the complete amplitude is squared 
to provide the event weight. 

The second approach is the CALKUL paradigm. In this approach, one 
concentrates on amplitudes with massless particles in the inital and final states, 
the typical situation for e^e~ annihilation when top quarks, W bosons, and 
other heavy particles have decayed to their final products. Then it is possible 



to compactly represent the amplitudfi for the whole process of production and 
decay in terms of spinor productsp\lB 

{12)^ul{pi)vr{p2) , m=UR{pi)vL{p2) . (1) 

The spinor products can in turn be computed directly for the set of initial and 
final massless four- vectors in a given event. 

It is important to note that there are no important polarization effects 
associated with hadronization, except that the r decay depends strongly on r 
polarization. This effect should be |accounted explicitly by decaying r's through 
the simulation program TAUOLAO 

The use of helicity amplitudes to systematically describe-LC physics was 
pioneered by the JLC group, using the programs HEliASliS for automatic 
Feynman diagram computation, and BASES/SPRINGli^ to provide weight- 
1 events. The current version of their package is PHYSSIM in Table El 

The need to build up systematically the full complexity of LC reactions — 
including beamstrahlung, initial-state radiation, initial- and final-state po- 
larization effects, and hadronization — has led the authors of almost all the 
simulation programs to embrace object-oriented programming for their future 
versions. One relatively simple generator, pandora, already segregates the 
beam and e^e~ reaction information into separate C++ classes aihich interact 
through a simple interface. Further details can be found in ref. lij. 

4 Supersymmetry 

The next few sections of this report will discuss generators devoted to specific 
problems of LC physics. The first of these is the coherent representation of 
supersymmetry processes. 

The full set of processes in e^e~ annnihilation to two supersymmetric 
particles is now available in three different prograrps., the supersymmetric ex- 
tension of PYTHIAP the lateslJ-elease of ISAJEtM and the SUSYGEN pro- 
gram written for LEP 2 studiesllH ISAJET was the first to include polarization- 
dependent cross sections. The new version also correctly includes the matrix 
elements for 3-body decays. SUSYGEN gives a complete treatment of initial 
and final polarization effects using the helicity-amplitude paradigm and even 
allows for nonzero phases in the A and gaugino mass parameters. ISAJET and 
SUSYGEN explicitly include parametrizations of beamstrahlung. All three 
programs allow input of a general set of supersymmetry parameters. Given 
the model-independent character of LC measurements, this is an important 
feature. PYTHIA and ISAJET also include facilities that compute the super- 
symmetry spectrum from an underlying modeK SUSYGEN includes an inter- 
face to spectrum calculations with SUSPECT ta 



5 Precision Standard Model 

For calculation of Standard Model background processes at a LC, it is not 
sufficient to consider e+e~ annihilation to on-shell 2-body final states. Back- 
grounds to new physics typically come from higher-order corrections in which 
additional fermions are produced or from e^e~ —>■ W~^W~ processes in which 
one W boson fluctuations far off the mass shell. In fact, these reactions are 
not distinct and one must include all e~^e~ —>■ 4 fermion Feynman diagrams in 
order to obtain a gauge- invariant result. 

This is already an issue at LEP 2 and a very serious effort has been made 
to provide 4-fermion event generators for the LEP 2 experiments. The status 
of generatorsJpr 4-fermion and W pair physics has recently been described by 
Bardin, et aZEj These programs typically use the CALKUL paradigm. Though 
their implemetations are slightly different, they agree excellently among them- 
selves and with the LEP 2 for configurations of four fermions all at large 
relative momepta. For breAtity, I have included_only three of these programs, 
EXCALIBURpa KORALWpi, and WPHACTp in Table |[ 

Two unresolved conceptual problems in the simulation of 4-fermion pro- 
cesses are the treatment of the W width for an off-shell W and the correct 
inclusion of transverse momentum for almost-coUinear initial state radiation. 
In both cases, there is no simple prescription which is gauge-invariant. This 
leads to discrepancies among the various programs in certain specific kinematic 
regions. For example, for the process e+e~ — > e^vdu (very low mass single 
W* production), the various generators give 10% differences in the predicted 
cross section when m{du) is as small as a few GeV. 

Additional challenges can be found in higher-order processes. For the study 
of the standard process e+e~ -^ it and also for many searches, one needs aii 
event generator for e^e~ — > 6 fermions. Accomando, Ballestrero, and PizziccJ 
have computed the relevant cross sections in a suitable form and are preparing 
a new generator, SIXPHACT. Alternatively, methods are now available to 
allow one to perform the computation automatically; I will describe these in 
the next section. 

At the same time, it is necessary to reach for higher accuracy in the siiikj 
ulation of 2-fermion final states. KORALZ, by Jadach, Ward, and WaspJ 
achieved an accuracy of 0.1% in the calculation of the small-angle Bhabha 
scattering cross section, and this was essential for the precision cross section 
normalization at LEP 1. At higher energies, it is necessary to treat multiple 
photon emissions coherently, and the precision calculations must be extended 
to larger forward angles. Jadach and collaborators have just released a new 
program KK which addresses these issues. 



6 Do it yourself 

Eventually, workers in LC physics will have a need for event generators for pro- 
cesses that have not been included in the standard programs. The traditional 
recourse in this case has been to find a friendly theorist with timejan his hands. 
Today, Japwever, amother course is made available by the GRAClES and COM- 
PHEPtJ programs. These encode the Feynman rules of the standard model 
and certain extensions, and allow encoding of arbitrary additional Lagrangian 
couplings, and then automatically generate the numerical sum of Feynman 
diagrams. The authors of GRACE have made available a supersymmetric 
extension which already includes all 239 e^e~ —^ 3-body supersymmetry pro- 
cesses, and all 1424 e+e~ -^ 3-body processes. A typical process might involve 
the summation of 100 tree diagrams. The calculations are done at the ampli- 
tude level, allowing the event generation to include full spin correlations using 
the helicity-amplitude paradigm. These programs can also compute one-loop 
corrections by evaluating diagrams in terma_^of the standard set of one-loop 
integrals defined by Passarino and VeltmanEII A more complete discussion pi 
these systems can be found in Perret-Gallix's contribution to this conferences 

7 Conclusions 

In this report, I have tried to summarize the array of programs that are now 
available to perform event generation for LC physics. These range from the 
general-purpose generators PYTHIA and HERWIG, to specific tools for super- 
symmetry and multi-fermion simulations, to tools for automatic generation of 
events for arbitrary physics processes. For the future, we expect to see trends 
toward object-oriented and modular programs, toward detailed high-accuracy 
computation of standard background processes, and toward further automa- 
tion of complex calculations. We are well on the way to the level of accuracy 
and generality that will be needed for the LC physics program. 
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